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Abstract 

Bethe-Salpeter equation is applied to p-p elastic scattering. The 
observables of spin are calculated in the framework of the M matrix 
using the two-body interaction potential. The parameter of the pseu¬ 
dovector coupling constant is adjusted so as to reproduce the spin 
singlet part. It is shown that the spin rotation R(0) and A{0) are 
improved by the resonance effect for 1 Sq. 

Elucidation of nucleon-nucleon two-body system in terms of the meson 
exchange model is the main purpose in the present study. The nuclear force 
between nucleons is mediated by the various mesons and their creation and 
annihilation are described by the method of the field theory. To relate the 
phenomenologically accepted ideas of the meson theory with the dynamical 
mechanism of two-body system the Bethe-Salpeter (BS) equation is indis¬ 
pensable and the investigation for a long time has revealed the applicability 
to various phenomena on the composite system. 

Recently our calculations for nucleon-nucleon elastic scattering have shown 
that the BS equation is suitable not only for bound state that is deuteron also 
for the scattering state [Tj. In the study the lowest-order calculation known 
as the Born approximation is corrected by the higher-order calculation and 
the results of the differential cross section and the polarization are improved 
to make possible to compare with the experimental data quantitatively. We 
proceed to calculations for various observables of spin in the elastic scattering 
as described later on. 

One of the characteristics in the present formulation by the BS equation 
is the inverse square potential. The virtual creation and annihilation process 
due to the mesons is represented by the Feynman propagator and for example 
the form of the a meson is expressed by using the modified Bessel function of 
the first order as ~ K\{r)/r in the coordinate space. Expanding the BS am¬ 
plitude by the set of the Gamma matrix and imposing the auxiliary condition 
the resulting equation appears to be equivalent to the Schrodinger eigenvalue 
equation with the potential terms as above on the various mesons. So, the 
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leading inverse square part (~ 1/r 2 ) is stronger than the screened Coulomb 
potential (Yukawa potential) at the short-range region and expected to have 
an effect on the two-body system. 

Another important character of the BS equation is that the BS amplitude 
consists of 16 expansion coefficients which can act as the wave functions of 
the relative coordinate between two nucleons. The role of each component is 
made clear by studying the property of spin. The set of three polar-tensor 
components is assigned the spin triplet (5=1) and therefore suitable for con¬ 
structing deuteron and the scattering state. Besides, the spin singlet (5=0) 
state is another object to be considered and also possible to understand it 
within the BS equation. It is the point at issue in the present formulation 
to select the most proper component among three, that is, the zeroth vector, 
the zeroth axial-vector and the pseudoscalar ones. While the zeroth vector 
component is independent of the other BS components in the simultaneous 
equations, the spin function /V of it is symmetric (*/V = TV) under the 
exchange between two fermions and then set identically equal to zero. On 
the remaining two components there is no reason which one is advantageous 
as the 5=0 wave and it is one of the main subjects to determine it by com¬ 
paring them with the experimental data. 

Difference between the zeroth axial-vector ( av ) and the pseudoscalar ( ps ) 
equations for the 5=0 part is represented by the strength of the leading order 
inverse square potential part Vi(r) = —g i M~ 1 /r 2 + 0(r°) (i = av or ps) given 
as 
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in which the isospin T = 0,1 of two nucleon system and M is the nucleon 
mass. For the isoscalar a and u mesons contribute dominantly, both g av and 
g ps are negative values and thus the interaction of the 5=0 part may become 
a repulsive force at the intermediate energy. The signs of g av and g ps are sig¬ 
nificant to construct the proton-proton (p-p) elastic scattering because the 
virtual (5, T, L) = (0,0,1) state is introduced to make the phase shift param¬ 
eter of the (1,1,1) state remains a real number [Tj. In the present calculation 
either g av or g ps dependence of the higher-order term is very large, then, to 
compare these two models we set g av = g ps and use the individual values of 
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the coupling constant g u of u meson. It is fortunate that the results of the 
calculation are not sensitive to g w because of the heavier value of the mass 
(mjm, = 5 ~ 6). 

In order to examine the meaning of each component of the BS amplitude 
the spin observables are appropriate to draw some conclusions and to find 
out the two-body interaction between nucleons. When one evaluates the spin 
observables it is formulated in terms of the density matrix[2] and the phe¬ 
nomena on the elastic scattering between two nucleons are represented by 
the quantities on the isospin T as 
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where the subscript /i of cd 1,2 ) for particle 1 or 2 is p, = 0, x, y, z and so on 

/-1 Q \ 

(do = 1). Extension of spin to the arbitrary directions is done similarly. 
Although the actual calculation of aJ 1VTp is restricted to which the number of 
the subscripts except 0 is at most two in the present study, the applications 
to the higher order quantities are possible. 

The elastic scattering is resolved by the method of the M matrix which has 
already been reported in our previous studyjl]. In it, to specify each partial 
wave of the M matrix the nuclear bar phase shift parametrization is used 
with no inclusion of the Coulomb force. The elements of the K matrix are 
associated with the standing wave solutions. In place of the exact solution 
the approximate form of the wave function under the leading order inverse 
square part of the potential derived from the BS equation is used in which 
the regularization procedure due to the cut-off function in the pion exchange 
interaction is taken into account. Our procedure to regularize the divergences 
caused by the inverse fourth power potential from the pseudovector coupling 
interaction acts only on the Born term, nevertheless, it has been verified that 
the results of calculation work well. 

A point to note is that in the present formulation the charge independence 
is broken. There are two reasons. For evaluation of the p-p clastic scatter¬ 
ing within the ladder approximation the pion-nucleon coupling constant / 
of the pseudovector interaction is reduced from the standard value of it to 
the ~ 1/20. Adding that the 3 P 2 state is corrected by the higher-order term 
using the inverse square part of the potential. Consequently the isospin T —1 
part of the M matrix (M 1 ) is different between p-n and p-p case through the 
phase shift parameters. 

Calculation of the M matrix is done in the center of mass system at the az- 
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imuthal angle p = 0 for the scattered nucleon. The elements of the M matrix 
are represented by the composite spin states |s>= 100> and |s z >= \ls z >. 
The irreducible representation of cd 1 ’ 2 ^ required is connected to the direct 
product representation by the unitary transformation. The matrix elements 
between spin singlet and the triplet Mj Sz and Mj zS vanish as the result of 
the parity conservation of the interaction for two identical fermions. The 
remaining 10 elements are reduced to Mj s , Mj 0 , and M^ 0 on 

account of the space reflection invariance of the M matrix. Consequently all 
spin observables of two nucleon elastic scattering are prepared as a function 
of the center of mass angle and the incident energy. 

The elements of the spin triplet are not independent and which are related 
to each other by the time reversal invariance relation[3]. To apply it in the 
present calculation the general relation of M{6 ip) is used at ip = tt. It is 
alternatively performed by the unitary transform on M{9 pi) to reverse the 
scattering plane before taking the trace with the term which breaks the time 
reversal invariance at p — 0. Then, it is verified that the condition of the 
time reversal invariance is satisfied in the numerical calculation. 

While our calculation with the BS equation is on the basis of the relativis¬ 
tic framework of the theory and the Gamma matrix expansion appears to 
maintain it, the transformation between the center of mass and the labora¬ 
tory system is done within the nonrelativistic kinematics for simplicity. It is 
conveniently performed by introducing the triad of the unit vectors ( P,K,n ) 
defined as P = (k + k')/\k+ k'\, K = (k' — k)/\k' — k\ and n = P x K — y, in 
which k (|| z) and k' are momenta of the incident and the scattered nucleons 
in the center of mass system respectively. The direction of the final state 
scattered nucleon (p s ) and that of the recoil nucleon (p r ) in the laboratory 
system is represented by the unit vectors as p s = P and p r = —K. The side 
directions of the final state scattered and recoil nucleons in the laboratory 
system are given by n x p s = K and n x p r = P accordingly. 

Investigating the difference between the av and ps equations for the ^=0 
part and the role on p-p elastic scattering the spin-correlation parameter 
C nn = ay y00 /I 0 (J 0 = ag 000 ) is useful because at the center of mass scattering 
angle 6 = 90° the quantity J 0 (l — C nn ) = \M] S \ 2 /2 is independent of the S =1 
of the M matrix and makes determine the singlet part M] s {6 = 90°). 

The quantity Jo(l — C nn ) of the p-p elastic scattering at 6 = 90° by the 
calculation is plotted in Fig. 1. The laboratory energy of the incident proton 
is 310 MeV. The experimental value of C nn is given by ~ 0.6 at the 300 ~ 400 
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MeV laboratory energy[1] and the related experimental value of / 0 (1 — C nn ) 
is 3.7x (1—0.6) ~ 1.48 which is shown by the dotted line. It is seen that the 
result of the av equation reproduces the M] s {9 = 90°) part at the intermedi¬ 
ate energy. The ps equation is about threefold as large as the experimental 
value and is not likely to reach the experimental area. 

It is necessary to examine the / dependence of the 5=1 part of the M 
matrix too. For the c o meson interaction is practically zero in the interaction 
of the polar-tensor equation giving the 5=1 wave, there is no model depen¬ 
dence between the av and the ps models. The quantity IqCrp = o x KPm 
calculated at 9 = 90° is plotted in Fig. 2. The spin correlation parameter 
Crp by the measurement is seen in the ref. [Sj. In order to adapt it to the 
definition of Crp here the relation Crp = —C_rp is used. The form of 
IqCrp = (| A/q x | 2 — |M 1 1 0 | 2 )/2 at 9 = 90° indicates the higher-order correction 
to the phase shift of the 3 P 2 state may influence the numerical value much. 
Comparing with the experimental value at the 400 MeV laboratory energy 
( IqCrp = 3.7 x (—0.32) ~ —1.18) the optimum value of / is around 0.05, 
a little smaller than the value of the 5=0 case. Using another result of the 
measurement j5] the experimental IqCrp value doubles and the point of in¬ 
tersection approaches the expected region. 

To investigate the value of / furthermore we take up the differential cross 
section J 0 which contains all the elements of the M matrix and has a basic 
meaning for the calculation of various spin observables. In Fig. 3 the / de¬ 
pendence is shown from 0.05 to 0.08 at 9 = 90° and the 310 MeV laboratory 
energy. The av equation gives the tentative value of / at ~ 0.065 well. Due 
to the overestimate of the M] s {9 = 90°) value the result of the ps equation 
is larger than that of the av equation. 

Next our interest is the angular dependence of various spin observables 
by using the value of the parameter / = 0.065. The differential cross section 
of p-p elastic scattering at the 310 MeV laboratory energy is plotted in Fig. 
4. The angular dependence is not seen noticeable to the difference between 
the ps and av equations. It implies that in the singlet part M} s (9) the 5- 
wave may play a dominant role compared with the higher waves even at the 
intermediate energy. 

The polarization P{9 ) = a 3 000 //o is suitable to observe the properties 
besides 9 = 90° since experimentally the maximum is observed at around 
9 = 30° ~ 40° and the relation P{9 ) = — P(ti — 9) due to the Pauli principle 
restricts the area under consideration. In Fig. 5 the polarization as a function 
of 9 at the 400 MeV laboratory energy is shown for the av and ps equations. 
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The quantity / 0 P($) is independent of M] s {9) and the underestimate of P{9) 
is connected with the overestimate of Iq{9) in the case of the ps equation. 
On the other hand the av equation gives the sufficient maximum value in 
exchange for appearance of the dip at 9 ~ 40° in the differential cross sec¬ 
tion. The magnitude and the position of the maximum value is understood 
by the enhancement of the 3 P 2 state due to the higher-order correction, that 
is, P(9) Pi (cos#) • Pl(cos9) rs-/ sin2# giving the position a rather shifted to 
the right than that of the experimental data[6j[7]. 

The depolarization D{9) = a y0y0 /Io covers the whole scattering angles 
and many experimental data have been accumulated nowadays[6j [8]. In Fig. 
6 the depolarization as a function of 9 at the 400 MeV laboratory energy 
is shown for the av and ps equations. The experimental D{9) values over 
6 = 90° are plotted by using the relation D{9 ) = K{n — 9), in which the spin 
transfer K{9 ) is defined as K{9) = a y00y /Io- At 9 = 30° ~ 40° the result 
of D(9) overestimates the experimental data as well as the av case in the 
polarization. In the present formulation only the P-wave is corrected by the 
higher-order terms making the partial wave be strong particularly[I]. Then, 
the extension of the method to include the other partial waves may improve 
the results of the calculation. 

In addition to the effect of the higher-order terms beyond the Born ap¬ 
proximation, resonant properties are thought to be essential for p-p elastic 
scattering at the intermediate energy. Particularly observables on the triple¬ 
scattering experiments contain M]. s (9) in the numerator and the correction 
may change the angular distributions largely. The spin rotation parameters 
R{9 ) and A(9) are concerned with rotation of the spin to the side direction 
K on the scattered proton initially directed to x and z and which are given 
by a} Km /I$ with i = x and z respectively. The calculations are performed as 
well as the procedure for D{9 ) and compared with the data of the measure¬ 
ments below 9 = 90° ill- 

In Fig. 7 the parameter R{9 ) is shown as a function of the center of mass 
angle 9 from 9 = 0° to 90° at the 400 MeV laboratory energy. It is clearly 
seen that both av and ps equations do not give the results convincingly. 
Although the difference between the calculation and the experiment is not 
small, it leaves room for changing the situation by taking into account the 
effect of the resonance. 

As the incident proton energy approaches the region of the pion produc¬ 
tion threshold the p-p system begins to scatter elastically via the process 
p + p —> n + -\- p + n primarily in the spin singlet 5-wave. Thus, the phase 
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shift parameter 6l (L— 0) is corrected by the resonant term and which is 
represented by 
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in which E and Eq are the incident energy of proton and the resonant energy 
located in the intermediate energy respectively^. The intermediate quasi¬ 
stationary state (7r + + p + n) has the width r associated with the life time 
r as r ~ 1/r. At present the definite value of r is not required because 
E ~ E 0 and so the second term of Eq. (4) gives ±7t/2 which acts on the 
element of the S matrix So = exp(2ih 0 ) to reverse the sign. 

The resonant state is expected to make effects particularly at 9 ~ 90° be¬ 
cause only three elements of the M matrix M] s) and M \ 0 remain therein 
and IqR(6) is shown by 


Iom = jWKWMliWl (0 = 90°). (5) 

The quantity I 0 A(9) is analogous to I 0 R(9) and it is given by replacing M^ x 
with —M \ 0 in Eq. (5). 

Taking into account the resonance mentioned above the av equation de¬ 
noted by av( 7t/2) nearly reverses the sign at 9 = 90° in comparison of the 
calculation without the resonance effect in the R{9). It makes possible to 
explain the trend of the experimental data from 9 = 10° to 50° well. The 
other resonance, for example, the P-wave is hoped to explain the trend at 
9 > 50°. For the calculation does not include the electro-magnetic interaction 
the sharp edge around 9 = 10° arising from the interference is not attained 
to describe here. While the approximate relation R(9 , h 0 + 7r /2) ~ —R(9, 5o) 
is seen in the case of the av equation, the ps equation does not show the 
feature which could be attributed to the difference of the phase shift of the 
higher waves (L > 2) in the spin singlet part. 

The resonance is applicable to the other spin parameters too. In Fig. 8 
the spin rotation A(9) is shown as a function of the center of mass angle 9 
from 9 = 0° to 90° at the 400 MeV laboratory energy. At 9 < 20° the results 
of calculation show the trend to decrease along with the experimental data. 
Without the effect of the resonance the results do not follow the sudden ris¬ 
ing observed in the measurements. By including the resonance the resultant 
curves denoted by av(n/2) and ps(tt/2 ) turn to increase as far as intersecting 
the experimental data at 9 ~ 70°. The remaining differences over 9 ~ 70° 
are not large in comparison of the case of R(9). 
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The improvement of R(0) by the resonance effect makes us apply it also 
to D{9) for satisfying the relation R(9) = D{6) at 0 = 0° regardless of lack of 
the electro-magnetic interaction. It has been verified that the overestimate 
at 0 — 30° ~ 80° in the calculation of D{6) is remedied by including the 
resonance effect. 

In the present study proton-proton elastic scattering is calculated by us¬ 
ing the Bethe-Salpeter equation. The change of the pseudovector coupling 
constant of pion is inevitable to reproduce the experimental data. Conse¬ 
quently the structure of the M matrix is modified and in fact the charge 
independence of the two-nucleon system is broken. The spin parameters in 
the triple-experiment are understandable if we take account of the resonance 
caused by the n + + p + n quasi-stationary state. 
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Figure Captions 

Figure 1: /o(l — C nn ) as a function of the pseudovector coupling constant / at the 
center of mass scattering angle 9 = 90° and the laboratory energy of 310 MeV. 
The av and ps denote the axial-vector and the pseudoscalar components for spin 
singlet part respectively. The experimental value 1.48 is from C nn data of 330 
MeV in ref. [3j and Jo value at 312 MeV and 9 = 30.67° in ref. J6j. 

Figure 2: IqCkp as a function of the pseudovector coupling constant / at the 
center of mass scattering angle 9 = 90° and the laboratory energy of 310 MeV. 
The experimental value -1.18 is from Ckp data of 400 MeV in ref. [5j and Iq value 
at 312 MeV and 9 = 30.67° in ref. [6j. 

Figure 3: Iq as a function of the pseudovector coupling constant / at the cen¬ 
ter of mass scattering angle 9 = 90° and the laboratory energy of 310 MeV. The 
av and ps denote same as that in Fig. 1. The experimental value 3.7 is from the 
data at 312 MeV and 9 = 30.67° in ref. [6]. 

Figure 4: The differential cross section Iq as a function of the center of mass 
scattering angle 9 at the laboratory energy of 310 MeV. The av and ps denote 
same as that in Fig. 1. 

Figure 5: The polarization P{9 ) as a function of the center of mass scattering 
angle 9 at the laboratory energy of 400 MeV. The av and ps denote same as that 
in Fig. 1. The experimental data are from the measurements at 392 MeV labora¬ 
tory energy in ref. [6] and 396 MeV in ref. [7]. 

Figure 6: The depolarization D{9) as a function of the center of mass scatter¬ 
ing angle 9 at the laboratory energy of 400 MeV. The av and ps denote same as 
that in Fig. 1. The experimental data are from the measurements at 392 MeV 
laboratory energy in ref. [6] and 366 MeV and 398 MeV in ref. [8]. 

Figure 7: The spin rotation parameter R(9) as a function of the center of mass 
scattering angle 9 at the laboratory energy of 400 MeV. The av and ps denote 
same as that in Fig. 1. The av{ 7r/2) and ps{ 7t/2) denote the axial-vector and 
the pseudoscalar components for spin singlet part respectively with the resonance 
effect. The experimental data are from the measurements at 392 MeV laboratory 
energy in ref. [6] and D u q s q at 398 MeV in ref. [8j. 
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Figure 8: The spin rotation parameter A{6) as a function of the center of mass 
scattering angle 9 at the laboratory energy of 400 MeV. The av and ps denote 
same as that in Fig. 1. The av( 7r/2) and ps(tt/2) denote same as that in Fig. 7. 
The experimental data are from the measurements at 392 MeV laboratory energy 
in ref. [6] and D^oko at 398 MeV in ref. [8]. 
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